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We investigate brane world models in higher-derivative gravity theories where the gravitational 
Lagrangian is an arbitrary function of the Ricci scalar. Making use of the conformal equivalence of 
such gravity models and Einstein-Hilbert gravity with a scalar field, we deduce the main features of 
higher-derivative gravity brane worlds. We solve for a gravity model that has corrections quadratic 
in the Ricci scalar and show one can evade both fine-tuning and the need for a bulk cosmological 
constant. An analysis of tensor and scalar perturbations shows gravity is localized on the brane and 
we recover the Newtonian limit. 
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I. INTRODUCTION 



Brane world scenarios envisage our universe as a hypersurface embedded in a higher-dimensional bulk. Large 
extra dimensions will not be visible in collider experiments if standard model fields are confined to this hypersurface. 
However, one should not expect gravity to be confined in the same way. Randall-Sundrum (RS) models [1, 2] are 
simple self-consistent scenarios in which gravity is effectively localized on a codimension one brane and the observed 
four-dimensional Newtonian limit is reproduced. One also recovers standard Friedman cosmology on the brane for 
energies low compared to the brane tension. The key is the non-factorizable, warped metric of the bulk. 

Nevertheless, RS models have an undesirable and much discussed feature: the problem of fine-tuning. For the RS 
II set-up which features one brane with positive tension A and a bulk cosmological constant A, we require 

A = -i^A 2 . (1) 

This fine-tuning is also necessary in cosmological models to obtain the correct expansion behaviour of the observed 
universe [3]. One of our goals here is to find ways to evade fine-tuning in brane world models and without needing a 
bulk cosmological constant. 

Since brane world models (mis)appropriate ideas from string/M-theory, it is usually hoped that they will arise as 
effective models in some string theory limit. From a phenomenological perspective therefore, one is led to include 
other fields, for example moduli fields, or to consider modifications to Einstein-Hilbert gravity. Bulk scalar fields have 
been extensively considered because they can stabilize two-brane models and also give rise to cosmological solutions. 
Because the Gauss-Bonnet term is nontrivial in five dimensions, the consequences of adding this term to the bulk 
action have also been explored. 

Here we combine ideas from both these extensions of the usual RS model. To be precise, we consider a gravitational 
Lagrangian which is nonlinear in the Ricci scalar: 
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S = -^ I d n x^\g\f(R). (2) 



We call this the physical or higher-derivative (HD) frame. In general, such an action leads to equations of motion of 
fourth order and there is little hope to solve them. However, this action is conformally equivalent to Einstein-Hilbert 
gravity with a bulk scalar field. In this scalar field (SF) frame, it is possible to obtain solutions to the equations of 
motion and to transform the results back to the physical frame. 

The outline of this article is as follows: first, we establish the conformal equivalence of the HD and SF frames, and 
place brane worlds within this framework. We present the complete action including the appropriate Gibbons-Hawking 
boundary term for /(i?)-gravity. In section III, we construct brane world solutions in the SF frame for potentials 
derived from superpotentials and for simple quadratic potentials. Finally, we perturb about these background solutions 
in order to show that gravity remains localized on the brane. Our methodology draws heavily on our recent work on 
bulk scalar field solutions [4]. We also include an appendix which details a calculation in the HD frame. 
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II. CONFORMAL EQUIVALENCE 

Higher-derivative gravity theories that depend only on the Ricci scalar are conformally related to usual scalar-tensor 
theories of gravity [5]. Consider the conformal transformation of the metric in n spacetime dimensions, 



Sab > £ 9AB = e "-f/w m (3) 

where is the gravitational constant of the full spacetime. The scalar field <f> acts like a dilaton field. Under this 
transformation, (2) becomes the Einstein-Hilbert action plus the usual scalar field action: 



S = d n x 



.±- R+ }. g AB (j)A(j)B _ v{4i) 



(4) 



The potential V(<j>) contains the information about the original gravitational Lagrangian. In fact it is related to the 
Legendre transformation of f(R). If we introduce tjj = f'(R) as the conjugate variable to R, then 

where 

f L ($) = il>R-f{R) and ^ = e Kn ^ ^ . (6) 

From the theory of Legendre transformations, we know that /l(V0 is well defined and concave as long as f(R) is 
concave. Then it will be possible to perform the inverse Legendre transformation to go back to the HD frame. 

Brane worlds, by their very nature, require us to add surface terms to the action. We introduce one codimension 
one brane with tension A as in the original RS II scenario and neglect matter fields on the brane. Then the brane 
action is 

S b = -xj d"-Vvl7|. (7) 
Under the conformal transformation (3), this becomes 

St = -\j rf"-V^HC/(</»), (8) 

where U(<f>) = e ~ K W^. 

To obtain consistent conditions on the boundary, we have to include the equivalent of a Gibbons-Hawking term to 
the surface action in the HD frame (see also [6]), 

Sgh = \ I d^aVhif'mKh, (9) 

K n J 

where K is the trace of the extrinsic curvature and [K}± = K + — K_. (We define the unit normal to the surface 
n A normal to point into the + region.) Transformation of the boundary term to the SF frame yields the familiar 
Gibbons-Hawking term. 

The equations of motion in the SF frame are 

Gab = 4 [Mb - 9ab {\g CD ^ c 4>D - V(<f>))] , (10) 

v 2 <f> = -V (</>), (11) 

subject to the junction conditions 

[Kab] ± = ^K 2 n \U(<P)q AB: (12) 

n A [<f> A } ± = -\U'(<j>), (13) 



where q A B is the projected metric on the surface[15]. 
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A. Relationship between f(R) and V((f>) 

While it is straightforward to compute V(<p) given an f(R), the resulting expression is, in general, rather unwieldy. 
However, in the situations we consider, V = on the brane. If we then assume that for small R, f(R) ~ -R, i.e. there 
is no bulk cosmological constant, then as an application of the mean value theorem, R = on the brane as well. It 
follows that in the vicinity of the brane 

1/(0) = iM 2 (0-^) 2 , (14) 

where 

*=y/&Wm and M = • [/ Wi ( ) • (15) 



III. BRANE WORLD SOLUTIONS 



From now on, and for obvious reasons, we specialize to the case of n = 5 dimensions. We make all quantities 
dimensionless by introducing the characteristic length scale L — k 2 / 3 . Then <j> — ► L~ 3 / 2 <f>, A — > L~ 4 A, etc. Actually, 
with a different rescaling, it is possible to absorb A as well, but this will make any fine-tuning of the parameters less 
apparent. In particular, we do not want the potential to depend on the brane tension. 

In the SF frame, where we perform all our calculations, the metric ansatz is 

ds 2 = e 2X ( v \ v dx' l dx v - dy 2 , (16) 

with the brane considered fixed at y = 0. For simplicity, Z 2 -symmetry is assumed; we will write down solutions only 
for y > 0. We consider perturbations about this background in the next section. 
It follows that the metric in the HD frame is 

ds 2 = e 2MY \ liV dx»dx v - dY 2 , (17) 

where 

A(Y)=X(y)-^(y) and Y = J\-^/ 2 ^dy', (18) 

so that Y = gives the position of the brane. Without loss of generality, we let A(0) = 0, which implies X(0) = 

The equations for the brane-bulk system are 

-3X yy -6X 2 = y 2 v + V, (19) 

6X 2 = i0 2 -V, (20) 

4> yy + 4X y <j> y = V, (21) 

with junction conditions 

x v (p) = -±\U@), MO) = \W{4>), (22) 

where we evaluate quantities as y { 0. It follows from the form of U and the constraint equation (20) that, as promised, 
V$) - 0. 



A. Superpotential approach 

Supcrpotcntials, familiar ingredients of supergravity theories, can be usefully employed for purely mathematical 
reasons [7-9]. Given a superpotential W(<j>), the associated potential in five dimensions is 

V = \W' 2 - \W 2 . (23) 
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FIG. 1: The gravitational Lagrangian f(R) for a linear superpotential with a 
negative for all a. 



-IX. 



The effective cosmological constant is 



(Smuggled into the definition of W is an overall factor of sgn(y), which we need to ensure Z2-symmetry across the 
brane, but this does not affect the calculation of V.) The advantage of such scalar field potentials is that the second 
order equations become first order and the constraint equation is immediately satisfied. Unfortunately, it is not 
possible to avoid fine-tuning in this class of models. Additionally, the transformation back to the HD frame leads to 
an f(R) that requires extra effort to make physical sensible. Nevertheless, these models are useful because we can 
write down exact solutions and they illustrate general features of brane world models in higher-derivative gravity. It 
is also possible to find models in which the effective bulk cosmological constant, A = ^/(0), vanishes. 
The system of equations becomes 

X y = -\W, <j> v = \W, (24) 
W@) = Ae" 2 ^, W'@) = --2=Ae- 2 ^. (25) 

Now it is clear that W = W((f>; a) — where a stands for the parameters of the superpotential — is a function of one 
variable that must satisfy two constraints. Thus we are forced to conclude that A = A(a), i.e. fine-tuning cannot be 
avoided. 



1. Linear superpotentials 

The superpotential W((f>) = 2(a0 -I- (3) leads to a quadratic potential. After application of the junction conditions, 
we have 

^) = -§a 2 (^-0)(V^-V3), (26) 

where a = — -^Ae -2 ^^. Note that the potential depends on A as expected and, because it is not of the form (14), 
we see immediately that the effective bulk cosmological constant cannot be zero. 

From this potential, we can, in principle, compute fj, (if)), R = f' L (tj}) and finally f(R). However, the inverse 
Legendre transformation does not exist when f'l(ip) = 0. This occurs at two points ip\ and corresponding to R\ 
and i?2, where ipi < < V>2 and R2 < R < < R\. Because R = f' L , as a function of Y, R(Y) will remain in the 
interval (-R2, Ri)- 

There are two options to make the theory physically sensible. First, we could imagine min{i?i, I-R2I} > Rp = 1, 
so that the theory of gravity as a whole would be modified before reaching these points. This is not so attractive 
because already \R\ ~ Ri, |i?2|- The second option is to assume that there is a second brane in the system at y = y*, 
equivalently Y = Y*. Subject to appropriate boundary conditions, 

X V {V*) = iA*E/(0*), Mv*) = -!A*f/'(^), (27) 
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FIG. 2: The gravitational Lagrangian f(R) for a quadratic superpotential that leads to a vanishing effective cosmological 
constant. In this case a = ^A. 



where quantities are evaluated as y j y*, this position can be chosen so that [tp, tp*] or [VVj "0] are contained in (ipi, ip2)- 
The use of a so-called regulator brane to 'slice off' singularities in a solution is well recognized in brane world models 
with scalar fields. Here, we are slicing off an offending piece of the theory space. In this picture, the two branes and 
the higher-derivative theory of gravity are inseparable; this is phenomenologically in the spirit of string theory. 

The expression for f(R) in the invertible range is not very illuminating; it is sketched in fig.(l) for a = — ^A. 
However, it is not hard to show that f(0) is negative for all allowed values of a. Therefore, as in the Randall-Sundrum 
case, the bulk cosmological constant is negative. 

The solution in the SF frame is straightforward to obtain: 



Hy) = <S> + a v, x (y) = jt^id + a v) - \ a2 y 2 , 



(28) 



so that 



A(Y) = -2 In 



1 + |Ae-^* /2 Y 



(29) 



Note that, in contrast to the original Randall-Sundrum case, A(Y) is smooth across the brane. This is a corollary of 
the junction conditions in the superpotential case and the form of U, but it is also a generic feature of higher-derivative 
theories. Because the equations of motion in the HD frame are fourth order, in general we expect discontinuities in 
third order quantities. It is sometimes thought that the jump in the warp factor is necessary to localize gravity on 
the brane, but we will show this is not the case in the next section. 



2. Quadratic superpotentials 

The extra degree of freedom of a quadratic superpotential allows us to construct models in which the effective 
cosmological constant vanishes. Letting W(<fi) = 2(a0 2 -I- /3(f) + 7) and adding the requirement that /(0) = to the 
junction conditions, we find 

n0) = -fa 2 (^) 2 (V-0-V3) 2 , (30) 

where now a — ^Xe^ 2 ^/^ . Near cj) — <p, we have V — — 2a 2 (<ft — <fi) 2 , which agrees with (14). 

As before, we have to be careful with the inverse Legendre transformation. The situation is as with the linear 
superpotential case except that R — 0. (Actually, f£ has additional zeros but there is only one interval (ipi,if>2) 
within in which we obtain a self-consistent solution.) Once more the final expressions are a little unwieldy and we 
resort to graphical representation: fig. (2) shows f(R) when a = ^A. There are +R 2 corrections to f(R) near R = 0. 

The solutions for <f>{y) and X(y) in the SF frame can again be straightforwardly written down and we omit the 
details. This time, however, there is no closed form expression for A(Y) in the HD frame. We depict it graphically in 
fig. (3), where it can be seen to be very smooth across the brane. 
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FIG. 3: The warp factor e A ^ Y ^ for a quadratic superpotential where A = and a = iA. 



B. Simple quadratic potentials 

The superpotential models are easy to solve and illustrative. In particular, we have learned that if f(R) has a term 
quadratic in R with positive coefficient, it appears we can avoid having a bulk cosmological constant. To show this is 
indeed a general result, we consider the simple quadratic potential, 

V{4>) = ±M 2 2 . (31) 

Using nonlinear perturbation theory, we recently [4] found brane world solutions for massive scalar fields in the 
limit of small M 2 and large — M 2 . (In terms of our scaling here, 'small' and 'large' mean with respect to A 2 .) Such 
potentials do not require fine-tuning, but the price we pay is in singularities that now crop up a finite distance from 
the brane. As indicated earlier, this generic feature of bulk scalar field models is remedied by the inclusion of a 
regulator brane. 

Returning to the HD frame is once again problematic since fi, is not concave everywhere — though where it breaks 
down is independent of M. We deduce from eq.(15), however, that in the region of invertibility 

f(R) =R- ^Ar 2 R 2 + .... (32) 

We indicate in the appendix that away from the brane this theory is quite different to a higher-derivative theory in 
which f(R) has only R 2 corrections. 



Small M 2 



Following [4], we let e = M 2 be our perturbative parameter and introduce the strained coordinate x = B 1 y. The 



zeroth order solutions in x are Xq 



i 

2 V / 3 V 



iln(l — lAa;), indicating a singularity at x 



details of the first order solutions X\ and <fii , we find 



C x th 2 M 2 
A(Y)^~-^M 2 / dxf^ = - T |g- — [2u 2 \nu(lnu~l) + u 2 -I], 
v Jo vo.x A 



Omitting the 



(33) 



where u = 1 — I Ax. Additionally, we must have B = 1 — jM 2 /\ 2 in order to avoid secular terms at first order. The 



coordinate Y in the HD frame is approximately 

Y = 2X- 1 



3B 
4-B 



1 



,3/ 4 



)• 



(34) 



Therefore, near the brane, A(Y) = ^AM 2 Y 3 + . . .. This result is in perfect agreement with the HD frame calculation 
given in the appendix and shows A(Y) will have discontinuities at the brane in its third derivatives. The sign of M 2 
determines whether the warp factor increases or decreases away from the brane. Since we want the latter case, we 
take M 2 negative from now on. 
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The situation as we approach the singularity is quite remarkable. While X and <fr both diverge, A is well behaved. 
To first order, and this is also borne out by numerical simulations, as u — > 0, Y — > Y s ~ 2A _1 and we have 

M 2 

A (X-) = m--tf <0 and A '( y «) = °- ( 35 ) 

Thus, while our solution in the SF frame fails at the singularity, it appears the solution in the HD frame will continue 
to Y > Y s . Although we have been unable to perform the HD frame calculation to prove this, it is moot since we 
require a regulator brane at Y+ < Y s to make the higher-derivative theory sensible. 



2. Large —M 2 

For large negative M 2 , the appropriate perturbative parameter is e = \M\~ . It turns out X and <fi functionally 
depend on both y and |Af|y. This means derivatives of these quantities are of lower order in perturbation theory. To 
first order, we find 

X(t/) = iln(l-iAy), <l>(y) = -^X\M\- 1 (l-^Xy)- 1 S m(\M\y). (36) 

There is a singularity as y — » 3A _1 but 4> is first order so remains small and the higher-derivative theory sensible 
until quite close to the singularity. At second order, X also picks up an oscillatory term, but A has an oscillatory 
component already at first order: 

A(Y) = | In (1 - | AF) + ^AlA/r 1 (l - ^F) -1 sin(|M|F). (37) 

This result is entirely consistent with the calculation in the HD frame outlined in the appendix. The smoothness of the 
warp factor across the brane and its oscillatory behaviour in the bulk is illustrated in fig. (4) for the case M 2 = — 100A. 

As pointed out in [4], the large negative M 2 limit in the SF frame looks like the usual RS model with vanishing 
bulk cosmological constant i/eq.(20) were ignored. It is this equation which gives rise to the fine-tuning condition 
(1). But the situation in the HD frame is even more curious: eq.(32) suggests f(R) — > R- In other words, in what 
looks to be Einstein- Hilbert gravity, we seem to be able to embed a brane with non-zero tension into a bulk without 
a cosmological constant! Of course, appearances are deceptive. For a start, this solution has undesirable limiting 
behaviour: while A — > Aq = | ln(l — gAF), Ay Aq,y etc. Second, f(R) — » R only in a certain range of R; we need 
a further brane to make the theory sensible. 



IV. PERTURBATIONS AND THE NEWTONIAN LIMIT ON THE BRANE 



The background solutions we have found for the higher-derivative theories do not have a kink at the brane in the 
HD frame. While it is known that a kink is not a prerequisite for gravity to be localized on the brane [10], i.e. in 
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order to have a zero mode, it is important to check that the brane-bound observer sees Newtonian gravity. Building 
on [4], we show that this is indeed the case in the framework of the simple quadratic potentials presented in the last 
section. 

For computational ease, we now choose a different gauge for the SF frame background: 

ds 2 = a 2 (y)(r] l , v dx»dx v - dy 2 ). (38) 
Scalar and tensor perturbations about this background are described by [11] 



SgAB = a 2 (y) 



(39) 



2A 



where is transverse and traceless, and indices are raised (lowered) by rf" {Vnv)- We have neglected vector 
perturbations because they will not have support to linear order. 

The tensor perturbations are already gauge invariant, and the pertinent gauge invariant scalar combinations are 

A = A-(B-E V ) V -^(B-E V ), (40) 
* = ^ + ^-(B-E y ). (41) 

Transforming back to the HD frame [11], we find anD = le - ^ 2 ^, Ahd = A + 7^tj<5></>, and ipHD = ip — j^^4>- 
However, B, E and are invariant. The gauge invariant quantities in the HD frame are 

Ahd = A + ^4> (gi \ (42) 
*hd = ^-^cj> (gl \ (43) 

where 5(f>^ = 5<f> + (f> y (B — E y ) is the gauge invariant perturbation of (j). Equivalently, we may write 

8^i) = ^ • [In f\R)\ SR i9l) . (44) 

We adopt the brane world equivalent of the longitudinal gauge in which B = E = from now on. 

In addition to the metric and scalar field perturbations, we must add a matter perturbation to the brane in the HD 
frame, £s m = ^■8m(l l HD)- The relation between the stress-energy tensors in the two frames is then 

r^=T(0)T« D , (45) 

where T(4>) = e^^l^ . Note that indices of the stress-energy tensor in the SF frame are raised (lowered) by ^ v (7^), 
where 7^^ = a 2 ^^. It follows that r = U(4>)thd- Actually, the distinction will not be important for us because the 
(j) appearing is the background <p on the brane and <j)(0) = <fi = in our scenario. Therefore a(0) = cihd(0) = 1- 

The final perturbation that must be considered is the perturbation in the brane position. We suppose the brane is 
perturbed from y — to y = C( a;A1 )- This scalar degree of freedom is manifestly gauge invariant in the longitudinal 
gauge and leads to an additional term in the perturbed junction conditions. We assume, however, that there is no 
matter on the regulator brane at y = y+, nor is its position perturbed. 



A. Tensor perturbations 

The tensor modes are the most important contribution to the Newtonian limit on the brane. It was shown in [4] 
that the Green's function for has the required 1/r behaviour and so we will not labour the point here. 

We found that there is always a zero mode of the perturbations h^ v which is proportional to a 3 / 2 , and that there 
are no tachyonic modes. Furthermore, for small M 2 and for large negative M 2 , massive modes of only make small 
corrections to the Green's function. The brane bending is given by □£ = — |r, as first pointed out in [12]. 
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B. Scalar perturbations 

We have to give special attention to scalar perturbations because they affect the transformation back to the physical 
frame. The danger is that the Newtonian potential may pick up an additional scalar mode in the HD frame. Our 
primary concern is the existence of zero and tachyonic modes. These have been carefully studied by Kofman and 
Mukohyama [13], and Lesgourgues and Sorbo [14]. 

The off-diagonal \p,, ^-components of the scalar equations of motion imply 

.4 = 2*. (46) 

The [5, 0]-component is a constraint equation, 

Wy + GH^ = -<p y 5(f> (9l \ (47) 



where H = a y /a, which we use to find 8(j)^ 9t ' . The one independent scalar equation of motion can be written as a 
wave equation: 



V yy - □* + - * y + 4 \H V - H^j * = 0. (48) 

The junction conditions at y = are 

H y 9i) (0) = ^Aa {u"8<t> {9i) - 2CJ'tf) + ^ aUr, „(0) = -±Aa {u'8<t> {ai) - 2tftf) . (49) 

The junction conditions at y = y+ are the same except for an overall sign and they do not include the term in r. 
This term comes from computing 8C$ m /8(j) in the SF frame, and it prevents the scalar perturbations from vanishing 
completely It also means we cannot fully adopt the machinery of [14]. 

We solve equation (48) by introducing a nonsingular Mukhanov variable. In our scenarios, H is never zero, so the 
variable that remains regular is v = a 3 ' 2 (84> {gi) - <f>y^/H). The equation for v is 



!ty -Hv-^fv = 0, (50) 



where z = a 3 ^ 2 (f) y /H. Since v y will have a jump at the brane positions, the equation for v can profitably be recast 
as a Green's function problem. Following the same procedure in [4], the problem reduces to finding the normalized 
modes ip m that satisfy 

(D + D_+m 2 )ij m = 0, (51) 
where D± = d y ± z v /z, subject to the boundary conditions d y ip m {0) = d y tp m (y*) = 0. It follows that 



m 2 > 



(52) 



Unfortunately, our ignorance of "0m at the branes prevents us from being able to compute the sign of the right hand 
side. We do know that z y j z becomes infinite if <p y — 0, thus the existence of tachyonic modes cannot be ruled out. 
However, the arguments of [14] suggest we can avoid them in the cases we consider. To be precise, if we neglect the 
matter perturbation on the brane, then 4> y not zero on [0, y*] implies [16] we have only modes with m 2 > 0. Wc now 
show explicitly that there are no zero modes. 
The general solution to eq.(51) for m 2 = is 

rv 

My) = ciz(y) + c 2 z(y) / Z - 2 {y')dy', (53) 
with ci, c 2 constants. If the boundary conditions are satisfied only if c\ = c 2 = 0, there is no zero mode. 

1. Small M 2 

In this case <j> y has no zeros. To a very good approximation, the background solutions give 

z(y)c2V3(l-^y /4 , (54) 
where r ~ 2 — |M 2 /A 2 to first order. It is straightforward to check that for all values of r, ip must be zero. 
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2. Large —M 2 



There is a zero almost immediately in (j> y , namely at y ~ ir/2\M\. To avoid this, the regulator brane must be 
brought very close to the visible brane, i.e. y+ ~ jMj" 1 . This is exactly the requirement we found in [4]. Then 



z{y)~2s/l(l-\M 2 y 2 ), (55) 



and once more ipo = 0. 



V. CONCLUSIONS 



We have initiated a study of brane world models in a class of higher-derivative gravity theories that are conformally 
equivalent to Einstein-Hilbert gravity with a scalar field. For illustrative purposes, we started with scalar field 
potentials that are derived from superpotentials. These exhibit the main features of higher-derivative gravity in 
brane world models. Then, we considered simple quadratic potentials which mimic the addition of R 2 terms to the 
gravitational Lagrangian. The typical background solutions are very smooth across the brane in the physical frame — 
in contrast to the usual scenarios. Furthermore, we found that we do not need a bulk cosmological constant and, 
therefore, can avoid the fine-tuning problem that plagues the original Randall-Sundrum scenario. Finally, we were 
able to ensure that gravity is effectively localized on the brane and that the Newtonian limit holds. Although the 
equations of motion are more difficult to solve in higher-derivative gravity theories, such theories offer an intriguing 
alternative to standard formulations of brane world models. We are now working on the cosmological consequences 
of this approach. 
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APPENDIX 



It is difficult to derive the equations of motion for f(R) theories of gravity, let alone to solve them. Here we present 
the simplest higher-derivative theory, 

f{R) = R + a 2 R 2 , (56) 

specialized to the case of Z 2 -symmetric brane worlds in five dimensions. 

The line-element in the HD frame that will give us both the equation of motion and the constraint equation when 
we vary with respect to the metric functions is 

ds 2 = a 2 {Y)r lllv dx^dx 1 ' - b 2 (Y)dY 2 . (57) 
Afterwards, we can put a(Y) = e A ^ and 6(F) = 1. In terms of dimensionless quantities, the total action is 

S = -i J d 4 xdYa%(R + a 2 R 2 ) + J d 4 xa 4 [(1 + 2a 2 R)[K]± - A] , (58) 

where 

4 (2a Y Y . 3a^ 2a Y b Y \ , r „, 8 ay 

R =Ta\ + ~t J—) and l K l± = -T ( 59 ) 

b l \ a a z ab ) b a 

The bulk equations are 

QA\r + 8a 2 {SAyyy A y - AA\ Y + 32A YY A Y + 5A Y ) = 0, (60) 

6A Y + 3A YY + 8a 2 (2A YYYY + 16A YYY A Y + 12A YY + 37 A YY A\ + 5A Y ) = 0. (61) 
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The junction conditions follow from taking the variations 6a, Say, Sayy etc to be independent on the brane. We find 

a 2 A Y (0) = 0, (62) 
6A Y (0) + 32(7 2 A Y yy(0) = -A, (63) 

and Ayy(0) = by Z2-symmetry. When a 2 ^ 0, Ay(0) — and so there is a jump in the third derivative only. 
Substituting (56) into (15), we find M 2 = - j|cr -2 and conclude A(Y) = -^XM 2 y' 3 in the vicinity of the brane. This 
confirms the result obtained in the SF frame. 

It is not possible to solve the fourth order equations in general. In the large a 2 limit, numerical simulations show 
that A decreases monotonically to a singularity at Y = Y s , where Y s is an increasing function of a. This is in 
stark contrast to the small M 2 limit away from the brane found in the text. On the other hand, for small a 2 , it 
is straightforward to show that (37) is indeed the approximate solution in the HD frame. We conclude that bulk 
singularities are likely to plague the higher-derivative theories considered here. 
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In the "mainly negative" signature, qab = Qab + tiatib and we define Kab = 5a C Vc«b and K — q AB Kab = V a ua- 
Technically, there are three additional criteria but they are met by our background solutions. 



